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A possible minimum toy model with negative differential capacitance for 
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We generalize a simple model for superlattices to include the effect of differential capacitance. It is 
shown that the model always has a stable steady-state solution (SSS) if all differential capacitances 
are positive. On the other hand, when negative differential capacitance is included, the model 
can have no stable SSS and be in a self-sustained current oscillation behavior. Therefore, we find 
a possible minimum toy model with both negative differential resistance and negative differential 
capacitance which can include the phenomena of both self-sustained current oscillation and I — V 
oscillation of stable SSSs. 
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I. INTRODUCTION 

Following the early pioneering studyi^ on vertical elec- 
tron transport in superlattices (SLs), one of the recent 
surprising discoveries is self-sustained current oscillations 
(SSCOs) under a dc bias 3,4 . A large number of exper- 
imental and theoretical studies have focused on differ- 
ent aspects of these oscillations. Experimentally, it is 
known that SSCOs can be induced by varying the dop- 
ing density^, temperature and magnetic field 4 ". Theo- 
retically, it is understood that SSCOs are accompanied 
by the motion of boundaries of electric field domains 
(EFDs) 5 . A discrete drift (DD) model capable of de- 
scribing both the formation of stationary EFDs and SS- 
COs emerged after many tedious analysises and numer- 
ical calculations^'?. Our understanding of SSCOs was 
greatly advanced through numerical investigations of this 
model 8 . However, this nonlinear dynamical model is too 
difficult to be solved analytically in a way that one can 
understand its different types of attractor solutions and 
their bifurcations. Thus, a minimum model with certain 
solvability may give some help for understanding these 
phenomena. 

Recently Wang and Niu(WN) 9 proposed a simple 
model with many nice features. This model is capa- 
ble of revealing analytically the source of instabilities in 
the electron transport. However, we can show that this 
model always has at least one stable steady-state solu- 
tion (SSS) under a dc bias. Thus it does not support a 
SSCO, and it cannot be a minimum model for the SSCO 
phenomenon. The question that we want to address is 
whether it is possible to find a minimum model which is 
simpler than the DD model. We shall show that a gener- 
alized WN model, by including the negative differential 
capacitance, may have no stable SSS, thus it can be a 
possible minimum toy model for SSCO in superlattices. 
This paper is organized in the following way. In SeclTTI 
we shall show that there is always a stable SSS in the 
original WN's model no matter what kind of I — V curve 
is used for each barrier. Then we will derive a generalized 



WN model by including negative differential capacitance 
in Sec lIIII In Sec llVI we shall show that this toy model 
with both negative differential capacitance and NDR can 
include both I — V oscillations of stable SSSs and SSCO 
behaviors if the model parameters are properly selected. 



II. EXISTENCE OF A STABLE STEADY-STATE 
SOLUTION IN THE ORIGINAL WN MODEL 

Let us first briefly review WN's model. WN employed 
the following three assumptions for their model of a su- 
perlattice of N quantum wells. (1) Inside each quantum 
well, the electronic states are described fully by time- 
independent quantum mechanics. (2) Charge carriers are 
in local equilibrium with each well so that a chemical po- 
tential can be defined locally. The chemical potential 
difference between two adjacent wells is call the bias V 
between them. (3) For a given bias V between two ad- 
jacent wells, a current I(V) passes through the barrier 
between them. In the following, Ij(t) and Vt(t) denote 
the electric current and bias through the i— th barrier, 
rii(t) the free charge in the i— th well. The current h(t) 
is supposed to depend on the bias Vi(t) only. WN thus 
obtain the following dynamical equations 



drii{t) 
dt 



= 7 i _ 1 (t)-7 i (t), 



where 



Vi(t)-Vi-i{t) = km(t), 



4-7T 

7s' 



(i) 



(2) 



(3) 



(4) 



Eq.Q and Eq.|J2J are the charge conservation law and 
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FIG. 1: The curve of I vs. V for barriers with a region of 

negative differential resistance in the case that both the min- 
imum current and maximum current of all barriers are the 
same. 

the Poisson equation for the i — th well, while Eq.|j3J) 
is the constraint condition that the total bias should be 
constant. As a first step, let us consider a simple case 
that I-V curves of all barriers are the same, as shown 
in Fig(T] Denote (I m ax,V max ) and {I m in,V min ) as the 
maximum and minimum currents and the corresponding 
biases in the I-V curve as shown in Fig^ An SSS of 
Eq.(|TJ— © 



should satisfy 



(5) 



(G) 



where / is the equilibrium current of the SSS. WN have 
proven the following conclusion 9 . 

Lemma 1. (a) An SSS is stable if all barriers at the 
SSS are in their PDR regions, i.e., 



dV 



Vl=V? 



Il(V?) > 0; 



(7) 



(b) if one and only one barrier of an SSS is in its NDR 
region, then the SSS is stable if and only if 



>0; 



(8) 



(c) if two and more barriers are in their NDR regions, 
then the SSS is unstable. 

With the use of Lemma 1, we can prove 
Claim 1. For Eq.QJ — J3J), there is at least one stable 
SSS for each value of U when all / — V curves are the 
same and has only one NDR region. 

Proof. When U G [0,NV max ] or U G [NV mm ,oo), an 
SSS with Vi = U /N for all barriers is stable according to 
Lemma 1(a) since V, = U/N drops in either lower- 
bias PDR region [0, V max ) or upper-bias PDR region 



(Vjnax, +oo). Thus our proof can be restricted to the re- 
gion of U G (NV max , NV m i n ), and all we need to do is to 
find a stable SSS for each value of U G (NV max , NVmin). 
Let us start by considering an SSS with the first barrier 
in its NDR region and all other barriers in their lower- 
bias PDR regions as shown in Fig^ The equilibrium 
current / of this SSS is slightly smaller than I m ax, i-e., 
51 = 1 — Imax < is a negative infinitely small value. 
The bias on the i— th barrier for this SSS is thus 



Vi(5I) = Vmax - SV 



with 



Then we have 



NV„. 



sv 

SI 



N 



= [I'M)}- 1 . 



(9) 



(10) 



N 



Y,v t (si) = siJ2&( v ^ 1 - ( n ) 



According to Lemma 1(b), this SSS is unstable when 
and only when 



A' 



which leads to 



N 



N 



NVrr. 



J2v t (si) = siJ2i I MT 1 >0- 



(12) 



(13) 



i=i 



Let us denote the total bias of the above SSS with the 
equilibrium current / as U(I = I m ax — SI) = 1 Vi (I). 
Now let us look for a stable SSS in the case when the 
total bias U = U(I). If the above SSS itself is stable, 
then it can be the SSS we look for. On the contrary, 
if the above SSS is unstable, then by Ea. (jl.Hf) we have 

U{I) = E l== i^( / ) < NVmax which means that U{I) 
is in the region [0, NV m ax\- Thus there is also a stable 
SSS in this case since we have shown that there is always 
one stable SSS when U € [0, NV max ]- Therefore, there 
is always a stable SSS when the total bias U = U(I) = 

J2iLi ^i(-O) regardless of whether the above SSS is stable 
or not. Since 51 is regarded as infinitely small, by Ea. i|ll|) 
the value U(I) — NVmax is also infinitely small. Thus we 
can say that we have proven the existence of a stable 
SSS when U S [0, U(I)]. Now let us further decrease I 
by SI and consider the corresponding SSS of the total 
bias U = U (I m cix — 251) with only the first barrier in the 
NDR region. Similar to Eqs.@-(|T3J), we can obtain that 
this SSS in unstable when and only when U(I m ax—2SI) < 
U(I m ax — SI). Therefore, we can prove the existence of a 
stable SSS when U G [0, U(I m ax — 2W)]. By decreasing I 
step by step, we can go on with this procedure before I = 
Imin ■ When / = Imin , the SSS we consider is as follows : 
the bias of the first barrier is V\ = V m in and the biases of 
all other barriers are Vi = Va- Thus we have proven the 
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existence of a stable SSS for U G [0, V A + (N - l)V min ]. 
Then, let us consider the following series of SSSs and 
increase the equilibrium current I from I m i n to I m ax '■ 
the first barrier in a part of its upper PDR region, i.e., 
Vi G (Vmin, Vb) and all other barriers in a part of their 
lower PDR regions, i.e., Vi G (VA,V max ). It is obvious 
that this series of SSSs are stable because all barriers are 
in PDR regions. Thus we have proven the existence of 
a stable SSS for U G [0, V B + (N - l)V max ] when I is 
increased to I m ax- At this point, the first barrier is in 
its upper-bias PDR region while all other barriers are at 
the cross point between the lower-bias PDR region and 
the NDR region. To go further, we can decrease I and 
keep the first barrier in its upper-bias PDR region, let 
the second barrier in its NDR region and all other N-2 
barriers in their lower-bias PDR regions. Following the 
above procedure, by decreasing / step by step we can 
arrive at an SSS that the biases of the first and second 
barrier are Vi = V2 = V m i n while all other N-2 barriers 
are at Vi = V A - Thus we have proven the existence of a 
stable SSS for U G [0, 2V A + (N — 2)V mm }. Then, follow 
the above procedure and consider the following series of 
SSSs and increase the equilibrium current I from I m % n 
to I max : Vi G {V min ,V B ) and V 2 G {V mmi V B ) while 
Vi G {VA,V max ) for all other barriers. Thus we have 
proven the existence of a stable SSS for U G [0, 2V B + 
[N — 2)V max ] when / is increased to I max - Then, we 
decrease /, and keep the first and second barriers in their 
upper-bias PDR regions, let the third barrier in its NDR 
region and all other N-3 barriers in their lower-bias PDR 
regions. Following the above procedure, we can prove the 
existence of s stable SSS for U G [0, 3V B + (N - 3)V max ] 
and so on so forth. Finally, we can prove the existence 
of a stable SSS for U G [0, NV max }. QED. 

Summarize the above proof, we can find the following 
general method to prove the existence of a stable SSS. 
Suppose that for some Uq > we have proven the con- 
clusion that there is always a stable SSS for U G [0, Uq]. 
Then, in order to extend the conclusion to a larger re- 
gion U G [0, Ui] with Ui > Uq, all we need to do is to find 
a series of SSSs that fulfill the following two conditions 




FIG. 2: The curve of I vs. V for (a) the first barrier, (b) 
the second barrier, and (c) the third barrier. The minimum 
current and maximum currents of them are not same. 



: (1) the maximum number of barriers in NDR regions 
for each SSS is one; (2) the total bias of the series of 
SSSs covers the region between Uq and U\. Using this 
method, we can extend the conclusion of Claim 1 to a 
general case and obtain 

Claim 2. For Eq.QJ — J3J), there is alway one stable 
SSS for each value of U when I — V curves of barriers are 
not necessarily the same while each barrier has only one 
NDR region. 

Proof. Since I — V curves of barriers are not nec- 
essarily the same, let us denote (I m ax,i,V max ,i) and 
(Imin,i,V m in,i) as the maximum and minimum currents 
and the corresponding biases in the I-V curve of the i— th 
barrier. Without loss of generality, let us suppose that 
Imax.i < Imax,2 < " < I m ax,N as shown in Fig|2I Let us 
start by considering a series of SSSs with the first barrier 
lying between the origin and C\ while all other barriers 
lying in their low-bias PDR regions. Let us increase the 
equilibrium current of these SSSs I continuously from 
zero to I m ax,i- This series of SSSs are stable because all 
barriers are in PDR regions. Now we arrive at an SSS 
with the first barrier at C\ while all other barriers in 
their low-bias PDR regions. Let us denote the total bias 
of this SSS as Uq. Then, we have proven that there is 
always a stable SSS for U G [0, Uo]- The rest task of the 
proof is to find a series of SSSs that fulfill : (1) the max- 
imum number of barriers in NDR regions for each SSS is 
one; (2) the total bias of the series of SSSs covers the re- 
gion (Uo, +00). Following the method in Claim 1, let us 
decrease / from I max ,i to I m in,i and consider a series of 
SSSs with the first barrier lying between C± and Di and 
the second lying between D2 and C2 while all other barri- 
ers lying in their low-bias PDR regions. Then we arrive at 
an SSS with the first barrier at D\ and the second barrier 
at C2 while all other barriers lying in their low-bias PDR 
regions. Let us denote the total bias of this SSS as U\. 
Thus we have proven that there is always a stable SSS for 
U G [0, U±]. Then, keep the first barrier in its high-bias 
PDR region, let / increase, and move the second barrier 
along the trace C2 — > D2 — > £2 —> F2 and other barri- 
ers correspondingly. If I m in,2 > Imin,i> then the second 
barrier can cross G2 when the first barrier is still in its 
high-bias PDR region. If I m in,2 < I7nin.i1 then the sec- 
ond barrier reaches F2 when the first barrier reaches D\ . 
At this point, the second barrier can not move continu- 
ously from F2 to H2 with the above method applicable. 
In order to fill the sudden jump from F2 to H2, we can do 
as follows : move the first barrier from D\ to B\ while 
all other barriers are unmoved. This decreases the to- 
tal bias. According to the above proof, there is a stable 
SSS for this total bias. Then, we can first move the sec- 
ond barrier from F2 to H2 while keep all other barriers 
in their low-bias PDR regions, then move the first bar- 
rier along B\ — > C\ — > D\ while keep the second barrier 
in its high-bias PDR region and other barriers in their 
low-bias PDR regions. By doing this, we let the second 
barrier cross its minimum current point while keeps the 
number of NDR barriers in this series of SSSs less than 
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or equal to one. 

Now let us consider the third barrier. If I m i n ,3 is larger 
than 2mm, i and I m in,2, then it is easy to let the third 
barrier cross its minimum current point. If 7 m in,3 is less 
than one of I min .i and I mm ,2, say,/ mm , 3 < J TO ,„,i while 
Imin.3 > Imin.2 as shown in Fig|21 then we can use the 
method in the proceeding paragraph to let the third bar- 
rier cross its minimum current point. If I m in.3 is less 
than both I m i n ,i and I m in,2i then we can do as follows. 
When the third barrier reaches E% , move the first barrier 
from D\ to B\. Then, move the second barrier from Gi 
to B2 when the third barrier reaches F3. After that, let 
the third barrier cross its minimum current point while 
keeping all other barriers in their low-bias PDR regions. 
Then let the first barrier cross its minimum current point 
while keep the third barrier in its high-bias PDR region 
and the second barrier in its low-bias PDR region. When 
this is finished, let the second barrier cross its minimum 
barrier with the use of the above method. Thus the first, 
second and third barriers cross their minimum current 
point. Keep this go on, we can cross the NDR regions 
of all N barriers. This finishes the proof because the 
conclusion is obvious when all NDR regions are crossed. 
QED. 

In the above proof, we restrict to the case that there 
is only one NDR region in I — V curve of each barrier. 
It is straightforward to generalize the above proof for 
general I — V curves that may have a lot of NDR regions 
separated by PDR regions. In fact, one can show that the 
above method works if the following condition is fulfilled : 
each barrier always lies in PDR region when the absolute 
value of its bias is above some upper-limit value. This 
condition is generally fulfilled. Therefore, we can come 
to the conclusion that there is always a stable SSS in 
WN's original model no matter what kind of I — V curve 
is used for each barrier. 



III. A GENERALIZED MODEL WITH 
DIFFERENTIAL CAPACITANCE 

In this section, we shall derive a generalized WN model 
and introduce negative differential capacitance. In the 
generalized model, we consider each barrier as a narrow 
homogeneous dielectric material. Then, instead of using 
the Poisson equation, we make use of the Gauss' law for 
dielectric materials 



[Di(t) - A-i(i)]S = 47mi(t) 



(14) 



where Di(t) the electric displacement in the i— th bar- 
rier and S the transverse area of the superlattice. We 
assume that Di(t) only depends on the average electric 
field Ei(t) — Vi(t)/Li, L. L the thickness of the i— th bar- 
rier. Eq.QJ still holds since it is the charge conservation 
law. Eliminate rii(t) by Eq.QJ and Ea. ll4|) . we have 

C i {V i )^-C i ^ 1 {V i -i) C ^=I i -i{V i - l )-I i {V i ) (15) 



and 



where 



i=l 



(It 



Ci{Vi) 



S dD, 



S dDi 



4tt dV, AttL; dE, 



(16) 



(17) 



One can see that CiiVi) serves as the differential capaci- 
tance of the i— th barrier because dDi/dEi can be viewed 
as the differential dielectric function of it. When the 
conventional relation dDi/dEi = e is substituted into 
Ea. (|17|) . one returns to WN's original model. Thus the 
above model is a generalized version of WN's model with 
bias-dependent differential capacitance Ci(Vi). 
Take a SSS of Eq.JTSJ) and fTfif 



which should satisfy 

I i -i(}t 1 )=I i (y°) 
and N 



(18) 



(19) 



(20) 



Follow the standard analysis of linear stability analysis 
as WN didS, we take a time-dependent perturbation of 
the following form 



(21) 



Substitute it into Ea. l|15|) and (|16fl and take the linear 
terms of 8Vi(t), we have 



G 2 °A + G$Ai = C%_ x \ + G?_iA-i 



iV 



where 



r o _ dh 



Vi=V° 



(22) 

(23) 
(24) 



are the differential capacitance and differential conduc- 



tance of the i— th barrier at Vi = V, 
Eq. (|22|) has the following form 

A 

A,, = 



(0) 



The solution of 



(25) 



where A is a non-zero constant. Substitute it back into 
Ea. (|22|) . we find that a non-zero solution should satisfy 



iV 



1 



C°A + G? 



N 
i=l 



i/c? 



= (26) 



The SSS is stable if and only if none of Xi satisfying 
Ea.p%|l is positive. 
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FIG. 3: Curves of I and D vs. V for two special cases which 
can lead to SSCO. (a) Barriers with negative differential ca- 
pacitance and with no NDR region; (b) barriers with two bias 
regions of negative differential capacitance and with an NDR 
region in between. 

IV. EFFECT OF DIFFERENTIAL 
CAPACITANCE 

In this section, we shall discuss the effect of differential 
capacitance on stability of the SSSs. We shall first con- 
sider the case when all the differential capacitances are 
positive, i.e., Cf > 0. As shown in Ea. (|26() in Sec lIIII an 
SSS is stable if and only if none of A^ satisfying F(Xi) = 
is positive. In the case when Cf > for all barriers, let 
us denote A, = -G°/G° and A x < A 2 < • ■ • < \ N . WN 
found a method to analyze the root of F(X) = Q£, and 
it is easy to show that their method works when all Cf 's 
have the same sign. Following WN's method^ we obtain 
: (1) F(X) = has no negative root when > for all 
barriers; (2) when Ai < and A2 > 0, -F'(A) = has one 
negative root if and only if F(X — 0) > 0. Since G 4 ° > 
for all barriers, the above conclusions lead to the same 
conclusions mentioned in Sec[H] : (1) an SSS is stable if 
all barriers have PDR, i.e., G° > 0; (2) if one and only 
one barrier has NDR, then the SSS is stable when 



JY 



N 



F(X = 0) = ^[G?]" 1 = ^Tim )]- 1 > 0. (27) 

z=l i=l 

Therefore, the same proof as in Sec[H] leads us to the 
conclusion that in the generalized WN model of Eq. I|15|) 
and there is always a stable SSS in the absence of 
negative differential capacitance. 

Now Let us see what can happen when negative dif- 
ferential capacitance is introduced. It should be noted 
that since at present we cannot have a reasonable way 
to expect negative differential capacitance in the origi- 
nal superlattice system, the model in this case can only 
be considered as a toy model. However, this toy model 
may be a minimum model for the phenomenon of self- 
sustained current oscillation. We shall show by a simple 
example that negative differential capacitance can lead 
to a non-stationary state. Suppose that all barriers have 
the same I — V curves with no NDR region, i.e., Gi > 0, 



while D — V curves are different and have a region of 
negative differential capacitance as shown in Fig^a). It 
is obvious that in this case the system has only one SSS 
at V® — U/N. Let us consider the case when V® = U /N 
drops into the negative differential capacitance region of 
every barrier, i.e., Cf < for all barriers. Then, with 
WN's method it is easy to check that this SSS is unsta- 
ble if two of \i = -G?/G? > are not equal. Thus 
the system has no stable SSS and must be in an SSCO 
state. This means that the generalized model with neg- 
ative differential capacitance and without NDR can lead 
to an SSCO state. However, the model with only neg- 
ative differential capacitance is not enough because the 
I — V oscillation behavior of stable SSSs observed in ex- 
periments can only be included by NDR. The SSCO be- 
haviors can also occur when both negative differential 
capacitance and NDR exist in the model. For example, 
consider the case shown in Fig|3{b) where each barrier 
has one NDR region between two negative differential 
capacitance regions. In this case, the number of SSSs is 
more than one. One can show that all SSSs are unstable 
when the average Vi — U/N drops in the region of both 
PDR and negative differential capacitance. Therefore, 
the generalized WN model with both NDR and negative 
differential capacitance may be considered as a possible 
minimum toy model for self-sustained current oscillation. 

In a recent study, Sanchez et. ali have pointed out 
that WN's original model will always tend to a stable 
stationary state if it exists. However, their argument 
seems not fit for our generalized model with negative dif- 
ferential capacitance. Following their argument, the total 
current of each barrier in our model is 



Itotal = Q(V^ + Ii(Vi) 

at 

which leads to a series of equations for Vi 
dt 

Put them into 



[CiiVi)}- 1 [Itotal - Ii(Vi)}. 



we have 



I1 



total 



(28) 



(29) 



(30) 



EJG(v^)]- 1 ' 

Put it back into Eq. H29|) . we have the equations for Vi 
dVi I liiVi) - IjiVj) 



(31) 



dt ^{v^ACM. 



(32) 



Unlike the case in the Appendix of Sanchez et. in 
the presence of negative differential capacitance G,* < 0, 
the total current Itotal in Eq. I|31(l can be time-dependent 
and the equations of V^'s in Eas. (|32|l are coupled together 
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through the factor [Cj (Vj)]^ 1 . Thus the existence of limit 
circles and other types of SSCO behaviors is possible in 
this toy model. Although no parameter in this toy model 
represents directly for some external parameters in the 
original SL system such as well doping, temperature and 
transverse magnetic field, the presence of negative differ- 
ential capacitance and NDR may be viewed as the con- 
sequence of selecting appropriate values of these external 
parameters in the original SL system. Possible detailed 
relation between this toy model and the original SL sys- 
tem needs further study and is not the aim of this paper. 

In summary, our analysis shows that a generalized 
WN's model without negative differential capacitance al- 
ways has a stable SSS and thus cannot have SSCO be- 



haviors. This means that negative differential capaci- 
tance is necessary for the existence of SSCO behavior. 
We show by a simple example that the same model with 
negative differential capacitance can have SSCO behav- 
ior. However, the model with only negative differential 
capacitance is not enough, since NDR is necessary for 
the I — V oscillation behavior of stable SSSs observed 
in experiments^. This means that the generalized WN 
model with both negative differential capacitance and 
NDR can serve as a possible minimum toy model for the 
phenomena of both self-sustained current oscillation and 
I — V oscillation behavior of stable SSSs. 

GX acknowledges the support of CNSF under grant 
No. I0347I0I and the Grant from Beijing Normal Uni- 
versity. 
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